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On  the  Use  of  the  Inverted  Abel  Integral 
for  Evaluating  Spectroscopic  Sources 


by 
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ABSTRACT 

A  numerical  method  for  evaluating  the  inverted  Abel  integral 
employing  cubic  spline  approximations  is  described  along  with  a 
modification  of  the  procedure  of  Cremers  and  Birkebak,  and  an  extension 
of  the  Barr  method.  The  accuracy  of  the  computations  is  evaluated 
at  several  noise  levels  and  with  varying  resolution  of  the  input  data. 
The  cubic  spline  method  is  found  to  be  useful  only  at  very  low  noise 
levels,  but  capable  of  providing  good  results  with  small  data  sets. 

The  Barr  method  is  computationally  the  simplest,  and  is  adequate  when 
large  data  sets  are  available.  For  noisy  data,  the  method  of  Cremers 
and  Birkebak  gave  the  best  results. 


BRIEF 

Three  methods  for  evaluating  the  inverted  Abel  integral 
are  presented  and  their  sensitivity  to  noise  in  the  input  data  and 
resolution  of  the  input  data  are  evaluated. 


f 


Introduction 

When  spectroscopic  observations  are  made  of  a  flame  or  plasma, 
the  quantity  observed  generally  represents  an  integration  of  the  con¬ 
tributions  to  the  signal  by  species  distributed  across  the  source  in 
the  line  of  observation.  Detailed  knowledge  of  the  conditions  at  any 
region  within  the  source  requires  the  application  of  special  tech¬ 
niques.  Two  line  atomic  fluorescence  may  be  used  to  obtain  local 
temperatures  within  flames  (1,  2).  The  droplet  injection  technique 
developed  by  Heiftje  and  Malmstadt  (3)  has  been  shown  to  be  useful 
for  making  spatially  resolved  temperature  measurements  (4).  Two  photon 
excitation  by  pulsed  lasers  has  been  proposed  by  Measures  as  an  instru¬ 
mental  technique  for  obtaining  spatial  resolution  (5).  These  tech¬ 
niques  require  special  instrumentation  and  are  generally  applicable 
only  to  the  study  of  temperature  profiles  or  distributions  of  species 
possessing  good  atomic  fluorescence  characteristics.  The  droplet 
injection  technique  may  be  more  general,  but  may  not  reflect  the  con¬ 
ditions  obtained  with  nebulizer  systems. 

Mathematical  techniques  may  be  used  to  extract  information  about 
the  interior  from  conventional  side  on  measurements.  When  applied  to 
sources  of  irregular  shape,  mathematical  techniques  require  that  measure¬ 
ments  be  made  from  more  than  one  direction  relative  to  the  orientation 
of  the  source  (6,  7,  8).  If  the  source  possesses  a  simple  and  regular 
geometry,  one  set  of  measurements  will  suffice.  This  is  a  common  case  in 
analytical  spectroscopy,  where  the  source  often  possesses  circular  symmetry, 
and  may  be  treated  by  the  inverted  Abel  integral  (9,  10,  11,  12,  13,  14). 
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This  integral  may  be  written  as: 


i  (r) 


1 

7 r 


rR 


r(x) 

-  r2)1'2 


dx 


(1) 


where  i(r)  is  the  value  of  the  variable  sought  as  a  function  of  the 
radius,  r;  I  (x)  is  the  value  of  the  observed  variable  as  a  function 
of  the  distance  between  the  optical  axis  and  the  source  axis,  x;  and 
R  is  the  outer  radius  of  the  source.  See  Figure  1. 

When  the  inverted  Abel  integral  is  used  to  determine  radial 
emission  profiles,  it  is  necessary  to  account  for  self  absorption  if  the 
flame  or  plasma  is  not  optically  thin  (15,  16,  17). 

In  this  laboratory,  we  have  been  interested  in  the  application  of 
the  inverted  Abel  integral  to  a  variety  of  sources,  including  inductively 
coupled  plasmas,  stable  analytical  flames,  and  transient  flames  produced 
by  pyrotechnic  flares. 

Because  the  transformation  of  a  laterally  observed  profile  into 
a  radial  profile  tends  to  accentuate  any  uncertainty  in  the  observations, 
it  is  important  to  evaluate  the  noise  sensitivity  of  any  technique  used 
to  perform  the  transformation.  An  extensive  analysis  of  error  propa¬ 
gation  in  Porter's  method  has  been  published  (18),  but  many  of  the  other 
methods  have  not  been  subjected  to  such  scrutiny. 

In  this  communication,  we  present  the  results  of  studies  into  the 
error  propagation  of  three  methods  for  evaluating  equation  (1),  The 
methods  investigated  were  a  modification  of  the  technique  of  Cremers 
and  Birkebak  (13),  an  extension  of  the  method  of  Barr  (12),  and  a  new 
method,  based  upon  Cubic  Spline  functions,  which  was  developed  by  the 
authors.  The  methods  are  compared  with  regard  to  their  accuracy 
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under  varied  conditions  of  data  resolution  and  noise  level,  and  recommenda¬ 
tions  are  developed  for  which  method  is  most  suitable  for  a  given  quality 
of  data. 

Description  of  the  Methods 

Application  of  the  inverted  Abel  integral  to  spectroscopic  data 
requires  the  construction  of  a  suitable  function,  I(x),  giving  the 
observed  quantity  in  terms  of  the  distance  between  the  source  axis  and 
optical  axis.  A  suitable  function  will  possess  the  following  properties: 
its  derivative  will  be  integrable  in  (1),  its  first  derivative  will  have 
a  value  of  zero  at  the  source  axis  ( i . e . ,  I ' ( 0 )  =  0),  and  it  will  reject 
noise  in  the  observed  data.  The  requirement  that  I ' ( 0 )  =0  is  dictated 
by  the  assumption  of  circular  symmetry,  and  a  failure  to  adhere  to  this 
constraint  will  generally  cause  a  large  error  in  the  computed  value  of 
i(0).  If  the  function  follows  noise  in  the  input  data,  the  resulting 
large  values  of  1 1 " ( x )  |  will  cause  significant  errors. 

Several  workers  have  used  polynomial  approximations  to  the  unknown 
function,  I ( x )  (11,  12,  13).  Such  functions  are  easily  generated  from  the 
observed  data,  and  their  derivatives  yield  an  integratable  representation 
of  Equation  (1 ) . 

To  allow  treatment  of  a  variety  of  curve  shapes,  the  data  are 
generally  broken  into  several  segments,  each  covering  a  portion  of 
the  range  of  x,  and  a  separate  polynomial  is  generated  for  each 
segment.  These  polynomials  are  spliced  together  to  obtain  I(x).  For 
the  jth  segment,  the  polynomial  representation  of  I ( x )  may  be  written: 

Ij(x)  ■  A0j  +  A,jX  ♦  A2jx2  ♦  43Jx3  *  ...  (2) 
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with  the  corresponding  derivative 

IJ'(x)-Au+2A2jx*3A3Jx2+...  (3) 

For  an  approxiamtion  of  this  form,  i(r)  is  obtained  by  a  summa¬ 
tion  of  Equation  (1)  over  the  appropriate  segments.  The  following 
equation  may  be  used  for  polynomials  of  fourth  or  lower  order: 


R  +  S  • 

V  Jr!'  +  SJt1  +  2fl2j  <Sj.l  -  Sj>  + 
j  j 


1  A 

2  rt3j 


2  ^i+1  +  ^’+1  1 

VlVl  -  Rjsj  +  r  *"  V  +  Sj  J  + 

4A4j  [S  J*V  S  J  +  ^  <Vl  -  sj)j 


(4) 


i.  l_ 

where  R.  is  the  inner  radius  of  the  j  data  segment,  or  r,  which- 

J 

ever  is  greater;  and  S.  =  (R2.  -  r2)1//2. 

J  J 


Development  of  the  Spline  Method 

Cubic  Spline  functions  consist  of  a  series  of  third  order  poly¬ 
nomials,  each  one  covering  an  interval  between  two  points.  These  poly¬ 
nomials  are  determined  so  that  the  composite  function  produced  by 
joining  them  end-to-end  has  the  following  properties:  the  function 
passes  through  all  of  the  points,  and  is  therefore  continuous;  the 
first  derivative  is  continuous;  and  the  rate  of  change  of  the  slope  is 
minimized,  that  is: 

/f" (x)2dx  (5) 

is  minimized  over  the  range  of  x. 

Cubic  Splines  have  been  used  for  interpolating  data  (18,19,20  ); 
however,  their  properties  of  having  a  simple,  continuous  first  deriva¬ 
tive  and  a  minimized  rate  of  change  of  slope  make  them  attractive  as 
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approximations  of  I(x).  Reinsch  (21)  has  addressed  the  problem  of 
fitting  Cubic  Splines  to  noisy  data  by  introducing  a  smoothing  parameter, 
S,  and  allowing  the  function  to  pass  outside  of  the  data  points,  subject 
to  the  condition 


n  (y j  -  fCx^)2 
i=l  (dyi )2 


(6) 


Here  y..  is  the  input  value,  f(x^)  is  the  value  computed  from 
the  spline  function,  and  dy.  is  the  weight  of  the  i L'  point.  N 
is  the  number  of  data  points.  Reinsch  suggests  using  the  standard 
deviations  of  the  data  points  as  the  weights,  dy^ ,  and  values  of  S  in 
the  range 

N  -  /ZTT  <.  S  _<  N  +  JZU 

Smoothing  values  well  in  excess  of  the  recommended  range  have 
been  explored,  and,  while  good  results  were  sometimes  obtained, 
performance  tended  to  be  erratic  when  smoothing  values  larger  than 
those  recommended  were  used. 

The  spline  method  is  implemented  by  first  normalizing  the  x 
values  to  an  outer  radius  of  1.0.  Reinsch's  algorithm  is  then  applied 
to  generate  a  smoothed  Cubic  Spline  function  having  N-l  segments.  A 
transformation  of  the  function  is  performed  to  give  it  coordinates 
compatible  with  Equation  (4).  Values  of  i(r)  are  computed  by  evaluating 
the  first  three  terms  of  Equation  (4)  over  each  segment  from  r  to  R. 

Oetails  of  the  transformation  are  given  below. 

Reinsch's  algorithm  produces  the  Cubic  Spline  in  terms  of  offsets, 
h.,  so  that  each  segment  starts  at  zero: 

J 

fj(x  -  hj)  =  A0j  ♦  A,j(x  -  hj)  +  A2j(x  -  h3)2  +  A3j(x  -  hj)3  (7) 


where  h.  is  the  value  of  x  at  the  inner  end  of  the  j  n  segment. 

J 

To  simplify  the  application  of  Equation  (4),  the  spline  equations 
are  converted  into  functions  of  x  as  follows: 


: 

A*0j  =  A0j  ■  Aljhj  +  A2jhj2  +  A3jhj3 

(8) 

A*lj  =  Alj  ‘  2A2 jhJ  +  3A3jhj2 

(9) 

A*2j  =  A2j  '  3A3jhj 

(10) 

A*3j  a  A3j 

01) 

This  yields  functions  of  the  form: 

f*j(x)  ♦  A*oj  +  A*,jX  +  A*2jx2  ♦  A*3jx3. 

(12) 

where  the  h j ' s  have  been  removed. 

One  further  modification  is  necessary,  since  the  condition 
I'(0)  =  0  is  not  met  by  the  Cubic  Spline  function.  The  innermost 
segment's  coefficients  are  recomputed  so  that  1 1 ( 0 )  =0  while  the 
original  values  at  the  ends  of  the  segments  are  preserved,  along 
with  the  continuity  of  the  first  derivative.  This  is  accomplished  by 
Equations  (13-16): 


00 

=  A*00 

10 

=  0 

3  (A 

20 

30 

II 

O 

■y^ - ^ - A*-j  -|  (N-l )  +  A31  (15) 


30  v"01  "  00M"  ' '  "  30'  ' 

These  yield  the  inner  segment  equation: 

f  *(x)  3  +  A#  x2  +  A*  x3 

T0  '  ;  00  M  20*  30 
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Note  that  the  value  Aq^  in  Equation  (16)  is  the  coefficient  from 
Equation  (7),  while  all  other  A‘s  in  the  right  hand  sides  of  Equations 
(13-16)  are  from  Equation  (12). 

U 

The  transformed  polynomials  frTg(x)  and  f*^(x)  through  f*fJ_i (x) 
are  used  in  Equation  (4)  to  compute  the  values  of  i(r). 

It  was  found  that  this  method  produced  large  errors  in  the  values 
of  i(0)  due  to  a  sharp  change  in  curvature  in  the  inner  segment  of  the 
spline  function.  The  following  modification  was  found  to  reduce  the 
error,  and  was  incorporated  into  the  method:  The  point  next  to  the 
center,  i.e.  at  r  =  j^y,  is  reflected  across  the  center  to  r  =  j(^y. 

The  spline  curve  is  then  fitted  to  the  N+l  points  from  r=  ^y-  to 
r  =  1.0.  The  segment  of  this  spline  function  covering  r  =  ^y  to 

r  =  0.0  is  discarded,  and  the  procedure  is  continued  beginning  at 
Equation  (8).  The  inclusion  of  this  extra  point  appears  to  steady  the 
function  near  r  =  0.  Reflection  of  2  or  3  points  was  also  tried,  but 
showed  no  benefit  relative  to  reflecting  1  point. 

A  further  attempt  at  reducing  the  error  in  i(0)  was  made  by 
reducing  the  statistical  weight,  dy,  of  the  points  at  r  =  ^y,  r  =  0 
and  r  =  yyp  This  technique  did  not  improve  the  results,  and  was  not 
incorporated  into  the  method. 

The  Cubic  Spline  method,  with  one  reflected  point  and  equal  weight 
on  all  points,  will  be  referred  to  as  SPLINE. 

Extension  of  the  Barr  Method 

Barr  (12)  developed  a  method  of  performing  the  Abel  inversion 
which,  while  using  a  mul ti segmented  polynomial  approximation  to  I(x), 
differs  from  other  methods  in  that  the  integration  is  performed  before 


-8- 


the  differentiation.  This  reversal  of  the  order  of  operations  is 
claimed  to  reduce  the  sensitivity  to  noise  in  the  input  data.  In 
addition,  this  method  greatly  reduces  the  amount  of  computation 
required  through  the  use  of  a  table  of  coefficients  from  which  a  sum 
of  products  is  used  to  compute  i(r): 
k  =  (N-l  )r 


i(r) 


1 


N 

(N-l  )7rn=k-2 


^knV 


for  k  >  2 


i(r)  = 


<  2 


(18) 


Note  that  i(r)  may  be  obtained  only  for  values  of  r  which  yield  inte¬ 
gral  values  of  (N-l)r,  and  that  the  function  I(x)  need  not  be  generated, 

since  only  the  observed  values,  IQ  through  In,  are  needed. 

The  table  of  coefficients  is  given  in  Ref.  (12)  for  N  =  21 .  Smaller 
data  sets  are  handled  by  the  assumption  that  I(x)  =  0  for  x  >_  1.0. 

In  order  to  perform  a  complete  comparison  of  this  method  to  the 

others  studied,  it  was  necessary  to  extend  the  table  of  coefficients 
beyond  N  =  21 .  The  table  has  been  extended  to  cover  N  £  51  by  the  use 
of  the  equations  presented  in  the  appendix. 


Modification  of  the  Cremers  and  Birkebak  Method 

A  method  involving  mul tisegmented  polynomials  in  which  the  segments 
are  overlapped  to  provide  smoothing  was  developed  by  Cremers  and  Birkebak 
(13).  The  method  has  been  implemented  in  this  laboratory  with  two  changes 
In  order  to  comply  with  the  restriction  that  I ' ( 0 )  =  0,  Cremers 
and  Brikebak  changed  the  form  of  the  approximation  of  I ( x )  from  Equation 

o 

(2)  to  a  fit  in  x  for  the  innermost  segment: 
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Vx>  =  A00  +  A10x  +  A20x  +  A30x  +  A40x8 

To  avoid  changing  the  form  of  Equation  (4)  to  allow  for  the  high 
orders  involved  in  (19),  we  have  used  Equation  (2)  throughout,  subject 
to  the  limitation  that  =  0  in  the  innermost  segment.  This  change 
from  the  original  method  might  also  be  expected  to  improve  the  results 
by  reducing  the  ability  of  the  inner  segment  to  reproduce  noise. 

The  amount  of  overlap  between  adjacent  segments  utilized  in  fitting 
the  approximations  to  I(x),  fixed  at  three  points  in  Cremers  and  Birkebak's 
work,  was  made  variable  in  our  modification,  and  is  expressed  in  percent. 
Thus,  a  50  %  overlap  specifies  that  one  half  of  the  points  from  the 
adjacent  segments  are  to  be  used  when  fitting  I(x).  This  modification 
permits  some  variation  in  the  amount  of  smoothing  performed. 

It  was  found  that  double  precision  (64  bit)  math  was  necessary  in 
the  least  squares  curve  fitting  computation  when  51  point  data  sets 
were  used. 

This  modification  of  the  method  of  Cremers  and  Birkebak  will  be 
referred  to  as  C  &  B  3  when  third  order  approximations  to  I(x)  are 
used,  and  as  C  &  B  4  when  fourth  order  fits  are  employed.  Except  for 
a  limited  set  of  experiments  intended  for  direct  comparison  to  Ref.  (9), 
all  work  with  this  method  involved  a  50  overlap. 
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Procedures 

Three  functions  describing  different  shapes  of  i(r)  vs.  r  were 
used  to  generate  test  data: 

i (r)  =  1  -  3r2  +  2r3,  (20) 

the  bell-shaped  curve  from  Ref.  (13), 

3 

i(r)  =  0.75  +  32r  over  the  range  0  5  r  5  0.25  and 
i(r)  =  16/27(1  +  6r  -  15r2  +  8r3)  over  the  range 
0.25  5  r  <  1.0,  (21) 


the  off  axis  distribution  from  Ref.  (13),  and 


1(r)  =  0.79788e'8r  , 


(22) 


a  Gaussian  curve  with  u  =0  and  c  -  k. 

"True"  inverted  data  were  prepared  for  comparison  with  the  various 
inversion  methods  by  evaluation  of  these  functions.  "Exact"  input  data 
for  the  inversion  calculations  were  prepared  by  integration  of 


rR 


I(x)  =  2 


x 


i(r)r 

- dr  (23) 

(r2  -  x2)1/2 


using  (20),  (21)  or  (22)  as  i(r).  These  data  were  rounded  to  +  5  x  10'^, 

corresponding  to  a  residual  peak  to  peak  noise  level  of  approximately 

-5 

0.02  ",  and  were  assigned  weights  of  5  x  10  for  use  by  SPLINE.  Various 
amounts  of  noise  were  added  from  a  normal  distribution  by  a  pseudorandom 
number  generator.  The  noise  was  expressed  in  percent  as  a  peak  to 
peak  value  referred  to  the  largest  observation  in  the  exact  data.  The 
standard  deviation  of  noisy  data  for  use  by  SPLINE  was  considered  to  be 
1/4  of  the  peak  to  peak  noise.  Data  sets  having  11,  21,  31,  41  and  51 
points  each  were  prepared  and  noise  was  added  at  0.1,  0.25,  0.5,  1,  2, 

5,  and  10  %.  Several  replicate  data  sets  were  prepared  at  all  added 
noise  levels. 

To  provide  a  direct  comparison  with  the  results  of  Cremers  and 
Birkebak,  data  sets  were  prepared  by  rounding  the  exact  values  to 
±  0.005. 

The  accuracy  of  the  inversion  calculations  was  evaluated  by  computing 
the  standard  deviation  of  the  errors  between  the  computed  inversion  and 
the  true  inversion.  This  was  expressed  as  a  percentage  noise  in  the 
inverted  data  relative  to  the  largest  value  in  the  true  inversion,  again 
using  the  criterion  that  four  standard  deviations  define  the  peak  to 
peak  noise.  Two  factorial  design  experiments  were  performed,  and  the 
data  were  subjected  to  Analyses  of  Variance. 

All  computations  were  programmed  in  FORTRAN  and  performed  on  a  Data 
General  NOVA  computer. 

Results 

The  modified  Cremers  and  Birkebak  method  was  compared  to  the 
original  using  the  off  axis  distribution  of  Equation  (21).  Both 


exact  data  and  data  rounded  to  ±  0.005  were  used,  as  in  Ref.  (.13),  The 
results  are  presented  in  Table  I.  The  modified  method  performed  much 
better  than  the  original  with  11  point  data  sets.  Above  11  points,  the 
original  method  was  better  when  exact  data  were  used. 

The  first  factorial  experiment  used  all  three  distributions,  11 
through  41  point  data  sets,  and  from  1  through  10  percent  added  noise. 

C  &  B  3  was  the  best  overall,  followed  by  SPLINE  and  then  BARR,  with  all 
differences  significant  (p<  0.05).  The  three  distributions  produced 
significantly  (p<  0.05)  different  results,  with  the  normal  curve 
(Eq.  22)  being  the  easiest  to  reproduce,  followed  by  the  polynomial  bell 
(Eq.  20)  and  the  more  difficult  off  axis  curve  (Eq.  21). 

Some  trends  are  evident  in  the  data  from  the  first  experiment  (see 
Table  II).  The  most  striking  is  the  failure  of  BARR  at  11  points.  Plots 
of  the  inverted  output  from  all  three  curve  shapes,  of  which  Figure  2b 
is  representative,  show  an  overdamped  appearance  when  compared  to  the 
output  of  the  other  methods  (Figures  2a  and  c).  Both  SPLINE  and  BARR 
appear  underdamped  in  the  presence  of  noise  when  a  large  number  of 
points  are  used,  as  depicted  in  Figures  3a  and  b.  At  the  lower  noise 
levels  (1  to  2  %) ,  BARR  with  31  points  worked  well,  as  did  C  &  B  3 
with  31  or  41  points.  At  higher  noise  levels  (5  to  10  %) ,  C  &  B  3  was 
clearly  superior  to  both  SPLINE  and  BARR.  All  three  distributions  gave 
qualitatively  similar  results. 

A  second  factorial  experiment  was  undertaken  to  elucidate  the  behavior 
of  the  methods  at  lower  noise  levels.  This  experiment  used  only  one 
distribution,  Equation  (21),  and  covered  the  0.1  to  0.5  %  noise  range. 
Statistical  analysis  grouped  SPLINE  and  C  &  B  3  together,  with  BARR 
significantly  (p  <  0.05)  worse,  due  again  to  the  failure  of  BARR  at  11 


-13- 


points.  The  results  of  this  experment  are  shown  in  Table  III.  All 
data  collected  on  the  Equation  (21)  distribution  are  summarized  in  Table 
IV  and  Figures  4  through  8. 

These  figures  show  that  SPLINE  provided  the  best  results  when  the 
input  data  contained  less  than  0.25  %  noise,  and  that  C  &  8  3  was  best 
with  data  containing  over  0.25  %  noise.  BARR  was  very  poor  with  11  point 
data  sets,  becoming  better  as  the  resolution  of  the  data  increased.  With 
small  amounts  of  noise  and  41  or  51  data  points,  BARR  approached  C  &  B  3, 
but  both  were  less  precise  than  SPLINE.  As  the  noise  level  increased, 

BARR  provided  results  intermediate  to  the  other  two  methods. 

As  seen  in  Figures  4  and  7,  C  &  B  4  was  inferior  to  C  &  B  3. 

Conclusions 

The  choice  among  these  methods  of  performing  the  Abel  Inversion 
depends  upon  the  situation  at  hand.  If  a  noise  level  of  over  0.25  is 
present  in  the  data,  the  best  results  may  be  expected  from  the  modified 
Cremers  and  Birkebak  method,  and  a  large  number  of  data  points  should  be 
used.  Where  precise  data  are  available,  the  Cubic  Spline  method  should 
provide  the  best  results,  if  a  reliable  estimate  of  the  standard  deviation 
of  the  data  is  also  available.  A  small  set  will  suffice  in  this  case. 

The  Barr  method  may  be  chosen  when  simplicity  of  computation  is  important, 
and  a  large  number  of  data  points  are  available. 
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Table  I.  Comparison  of  the  Mod. .ied  Cremers  and  Birkebak  method  to  Cremers  and 

Birkebak's  results  for  the  curve  given  by  Eq.  21.  All  data  are  peak  to 
peak  noise  as  a  percent  of  the  largest  value  in  the  true  inversion. 

aValues  have  been  converted  to  percent  noise. 

bThree  point  overlaps  were  produced  by  the  following  overlap  percentages: 

11  points  -  100  %■,  21  points  -  75  31  points  -  50  %,  41  points  -  37.5  % 

and  51  points  -  30  ".  Note  that  the  maximum  overlap  available  with  11  point 


data  sets  is  only  2  points. 


Table  II.  Noise  in  the  computed  inversion  for  levels  of  noise  in  the  input 
data  of  up  to  10". 

Each  value  is  the  mean  of  2  replicates. 
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Table  III.  Percent  noise  in  the  computed  inversion  for  three  methods 
with  input  noise  of  up  to  0.5  %.  All  data  are  for  the  off 
axis  distribution  of  Eq.  21. 

aEach  value  is  the  mean  of  5  replicates. 


M 

run 


0. 

0. 


0.1 

0.25 

0.5 

0.89 

1.52 

2.59 

0.49 

1.10 

1  .86 

0.79 

1  .17 

1 .92 

0.55 

1  .43 

1  .50 

0.58 

1  .75 

2.87 

25.60 

25.61 

25.64 

5.46 

5.56 

5.66 

2.33 

2.53 

2.77 

1 .36 

1.54 

2.14 

0.93 

1 .46 

1  .94 

2.6 


2 

5 

5.1 

8.8 

4.1 

19.9 

6.0 

13.4 

1  .78 

j  2.2 

1  .57 

;  2.2 

1  .35 

1  2.1 

1  .05 

*  2.2 

\ - 

!  ! 

1.46 

1 .70 

1.95 

I 

2.76 

1 

5.1 

4.9 

18.3 

36.8 

2.82 

2.88 

3.04 

2.87 

3.4 

4.3 

13.8 

22.2 

Table  IV.  Summary  of  data  collected  for  the  distribution  given  by  Eq.  21. 


Comparison  of  some  computed  inversions  with  the  true  inversion  of  the  Equation  21 
distribution.  Inversions  computed  from  51  point  data  sets  having  5 %  noise  added. 
Solid  line  =  true  inversion;  broken  line  =  computed  inversion. 

3a.  Inversion  computed  by  SPLINE. 
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e  6.  Noise  in  the  computed  Inversion  as  a  function  of  the  noise  level  in  the  input  data  for  31  point  data  sets  from 
Equation  21. 


ure  7.  Noise  in  the  computed  inversion  as  a  function  of  the  noise  level  in  the  input  data  for  41  point  data  sets  from 
Equation  21. 
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Igure  8.  Noise  in  the  computed  inversion  as  a  function  of  the  noise  level  in  the  input  data  for  51  point  data  sets  from 
,  Equation  21. 
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Appendix  -  Computation  of  3^ 

The  table  of  coefficients,  S.  used  in  the  Barr  method  mav  be 

Kn 

computed  for  a  system  having  X  +  1  points  by  the  following  algorithm  (22). 

The  table  is  general,  and  may  be  used  for  any  smaller  data  set. 

Copies  of  the  table  for  N=51  are  available  from  the  authors. 

For  k:  =  <t>,  step  1  until  X 
do  begin 

\  =  44S0k6  -  704k4  +  5352k ^  +  57.6 

C0  =  (-448k5  -  1088k4  ♦  96k-3  +  819.2k2  +  426.4k  -  2.4)/I<k 
Cj  =  (-224k3  +  544k4  -  1920k-3  -  212.8k2  -  158.8k  +  58.4J/K,. 

C2  =  ClOSSk4  -  1212.8k2  -  72)/Kk 

C.  =  (224k5  +  544k4  +  1920k°  -  212.8k2  +  158.8k  +  38.4)/K. 

C4  =  (448k3  -  1088k4  -  96k3  +  819.2k2  -  426.4k  -  2.4)/*^ 


For  n:  =  0,  step  .1  until  X 
do  begin 

IF  k  <  2  or  n  -  k  >  2 
4 


3ki>  Cl’kn 
1=0 


ELSE  IF  n  -  k  >  -2 
n-k+2 

kn  =  2  Ciak-2+1fn 

i=0 


end 


end 
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